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a supersymmetric standard model with bilinear R-parity violation. In this 
scenario, the mixing matrix among the neutrinos is in a restrictive form. We 
find that only the small angle MSW solution is allowed for the solar neu- 
trino problem when the results of the CHOOZ experiment are combined with 



the mass squared difference and mixing angle suggested by the atmospheric 
neutrino data. Collider signals of this scenario are also discussed. 
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Neutrino masses and their mixing are an intriguing window to physics beyond the Stan- 
dard Model. Neutrino oscillations among massive neutrinos have been strongly suggested 
by the atmospheric neutrino data at Super-Kamiokande ]I| as well as by the solar neutrino 
problem [Q. 

Here we shall consider this issue in the framework of a supersymmetric (SUSY) extension 
of the standard model |§. Usually in the supersymmetric standard model, one imposes Im- 
parity conservation to avoid too fast proton decay. However, it is well known that the 
stability of the proton does not forbid all R-parity non-conservation, but allows R-parity 
violation such that either lepton number or baryon number is broken, because the nucleon 
decays normally require not only the baryon number violation but also the lepton number 
violation. 

Actually it is widely recognized that the R-parity violation with the lepton number 
non-conservation |4]|| is a promising approach to give neutrino masses and mixing without 
introducing right-handed neutrinos ||. Here we shall consider the case where the R-parity 
is broken only by bilinear terms |J . This is particularly interesting because it is compatible 
with the idea of grand unified theory (GUT) where quarks and leptons typically belong to 
one representation of a GUT group. 

An important feature of this mechanism of generating neutrino masses is that one may 
be able to see the R-parity violation in future collider experiments [^,§ . 

In this paper, we shall investigate the nature of the neutrino mass matrix in this scenario. 
We find that large angle solutions to the solar neutrino deficit are not allowed when we com- 
bine the CHOOZ result with the v T -v^ solution of the atmospheric neutrino problem. 



This singles out the small mixing MSW solution |TT| for the solar neutrino in this scenario. 
We will then show that the soft SUSY breaking masses should be finely adjusted to simul- 
taneously explain the atmospheric neutrino anomaly and the solar neutrino problem. We 
will also discuss implications of this specific pattern of the R-parity violation at collider 
experiments. 

Let us first describe the model we are considering. The particle contents of the model 



are the same as those of the minimal supersymmetric standard model (MSSM). R-parity 
breaking bilinear terms are assumed to exist in superpotential 

W = fiH D Hu + HiUHv + Y?LiH D EZ + YPQiHuD* + YgQiHuUf. (1) 

Here H D , H v are two Higgs doublets, L t a SU{2) L doublet lepton, Ef is a singlet lepton, 
Qi a doublet quark, Iff, D c { a singlet quark of up and down type, respectively. Suffices i,j 
stand for generations. In this basis, the Yukawa coupling matrices for the charged leptons 
and for the down-type quarks are assumed to be diagonal with respect to the generation 
indices. The soft SUSY breaking terms in the scalar potential are taken to be of the form 

U soft = BH D Hu + BiLiHv + m 2 Ho H D H ] D + m\H v El + m^Z^ + ml^l) + ■■■, (2) 

where we have written only bilinear terms explicitly. In this paper, we assume lepton-Higgs 
universality in the soft masses, namely 

B % jB = ft/ ft, m 2 Lij = dijinl., m\. = m 2 Ho , m 2 HL . = (3) 

These conditions are achieved, e.g. in the minimal supergravity Q or in scenarios of gauge 
mediated supersymmetry breaking (GMSB). Our motivation to consider this universality 
is alignment of vacuum expectation values with the lepton number violating bilinear terms 



14| . We shall return to this point later. 

Of course, the lepton-Higgs universality is not stable against radiative corrections and 
we assume that Eq. (0) holds at an energy scale where these soft masses are given as the 
boundary conditions of the renormalization group equations. 

In this model, the R-parity violating terms are parameterized by 



1 Note that in the gravity-mediated-supersymmetry-breaking scheme, the lepton-Higgs universality 
in the soft terms may be violated through interactions relevant above the Grand-Unified scale p 
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Here, for simplicity, we have taken /x and /ij to be real. 

With (D) it is convenient to define a new basis of the Higgs and lepton supermultiplets 
in which the all bilinear mixing terms between the Higgs and leptons in the superpotential 
are absent: 
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We can neglect the Yukawa couplings for the first and second generations. This approxima- 
tion yields a global U(l) symmetry, resulting in one massless neutrino L 1 = CiL^ — S\L e . 

The rotation to take the new basis gives rise to R-parity violating trilinear couplings in 
the superpotential as well as the softly SUSY breaking terms in the scalar potential. For 
example a Yukawa coupling for bottom quark is written 



Y b Q 3 H D D c 3 = Y b Q 3 (c 3 H D - s 3 L 3 )D c 3 

Similarly a Yukawa coupling Y t L 3 Hdt c now reads 

Y T (c 2 s 3 H' D + C2C3L3 - s 2 L' 2 )(c 3 H' d - s 3 L' 3 )B c 
= s 2 s 3 Y t L' 2 L' 3 t c + Y t c 2 L' 3 H d t c - Y T s 2 c 3 L' 2 H' D r c 
~ Y T L T s 3 (s 2 (c 1 L fl + SiL e ))r c + Y t L t H d t c 
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Note that as far as S3 is small, the old basis in Eq. ([I]) approximates a mass eigen basis for 
charged leptons. On the other hand, the new basis with prime gives in good approximation 
a mass eigen basis for neutrinos, as we will see shortly. 



As we mentioned before, radiative corrections violate the relation (3), resulting in vacuum 
expectation values(VEVs) of sneutrinos. An analysis shows that v' z and u' 2 develop VEVs. 



In fact the renormalization group effects induce terms of the form [[HJ 

AB 3 L' 3 Hu + Am 2 HL3 L' 3 H% + AB 2 L' 2 H V + Am 2 HL2 L' 2 H%, (8) 

which gives the following sneutrino VEVs 

, M ABitcm8 + Am 2 HT , . % , 

M) = —i 1^2 ^^ ( H d) (i = 3,2) 
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Note that AB 3 and Am 2 HL3 are proportional to the bottom Yukawa coupling squared, while 
AB 2 and Am 2 HL2 are proportional to the tau Yukawa coupling squared, and thus the ratio 

tan0=)4f (10) 

is generally small. 

The non-vanishing VEV for the sneutrino induces mixing between the neutrino and the 
neutralinos, giving rise to a small mass for it due to the see-saw mechanism. Here it is useful 
to define a new basis 

v 'i = v 'z cos + ^2 srn 'Pi u 2 = u 2 cos — ^3 s i n 'Pi v '[ = v '\- (11) 

Then v 3 acquires a tree-level mass 

1 1 

mtree = t;9z v I XX^^W - s w O B ^f . (12) 
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Here m^o (i = l,---,4) represent neutralino masses and B ^o and O w ^o stand for the 

contamination in x° of Bino and Wino components, respectively. Also v v = J(v 3 ) 2 + (y' 2 ) 2 - 

The neutrinos acquire masses also at the one-loop level. The mass matrix induced 

from loop diagram contribution are diagonal for the v[ basis. Relevant diagrams we should 

consider here are those which include the trilinear R-parity violating couplings. A one-loop 

diagram including scalar bottom quarks induces a following mass for vL: 
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Here wf.jp denotes a left-right mixing mass squared in the sbottom sector, and m^ and 
m~ b2 are the eigenvalues of the scalar bottom quarks. Similarly a one-loop diagram involving 
scalar tau leptons gives a mass for i/ 2 
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Generically the latter contribution gives a dominant contribution to the mass of the second 
heaviest neutrino. The other neutrino remains massless in the approximation that the 
Yukawa couplings for the first two generations are set to be zero. 
Summarizing, the mass matrix is written in the v" basis 



m tree + m°, op cos 2 + m°° p sin 2 (m°, op — m°, op ) sin cos 
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When diagonalizing this matrix, the mixing angle for (i/ 3 , i/^plane is 
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(m tree + (m^ op - m^°, op ) cos 20) m tr . ee 

Because the tree-level mass is generically larger than the one-loop corrections, we find that 
the v'l basis is in a good approximation a mass eigenstate. 



Thus the mixing matrix of the neutrino sector, the MNS matrix [ll|, becomes 



U T 3 U r2 U T \ 



Us, 



\ 



U^ UfjQ u^ 

U e 3 U e2 U e i 

1 
c\ —s\ 
si ci i 



c 2 -s 2 
s 2 c 2 
1 



C0+<5 — s</,+<s 

S(/>+8 C^ + s 

1 



c e -s e 
-CiSfl cic e —S\ 
sis e sic e C\ 



(17) 



with 6 = 62 + <p + 5. Here i and a denote the weak current basis and the mass eigen basis, 
respectively. Note that the MNS matrix obtained is of a very restrictive form with only two 
angles. 

Now we would like to determine the mixing angles by using experimental data. Here we 
assume the hierarchical mass structure, namely m 3 ^> m 2 ^> m\ so that Am| 2 ~ Am^ ^> 
Amfi, as will be naturally realized in the model we are considering. 

Let us first consider the atmospheric neutrino anomaly, which is explained by an oscil- 
lation v^-Vr. The transition probability is 

P(u^u T )^4\U, 3 \ 2 \U r3 \ 2 sin 2 ^L 

,,222-2 ^ W 32 r /,n\ 

= 4c 1 s e c sin -j^L. (18) 



A recent analysis of the atmospheric neutrino at SuperKamiokande [17] indicates the fol- 
lowing ranges at 90% CL: 

Am 2 32 = (2 - 6) x l(T 3 eV 2 (19) 

Ac\s 2 e c 2 e > 0.85. (20) 

On the other hand, the CHOOZ experiment || gives a bound on the mixing angle 
U e3 = S1S0 to be 

4:sls 2 g (l - s\sl) < 0.2 (21) 

for the mass range suggested by the atmospheric neutrino problem. 

Combining Eq. (|20|) with Eq. (|2l|), one finds that the s± has to be small. This non-trivial 
relation among U e3 , U^ and U r3 in the MNS matrix restricts patterns of the neutrino oscil- 
lations. In fact, for small U e3 implied by small Si, the three flavor solar neutrino oscillation 
is well approximated by the usual two flavor oscillation solutions, namely the small angle 



MSW, the large angle MSW, and the vacuum oscillation fllQ |. However, the latter two are 
excluded because the mixing angle involves s\. This argument enables us to single out, in 
this model, the small angle MSW solution as the only possible explanation in terms of the 



neutrino oscillations to the solar neutrino problem. This is one of the main results of our 
paper. 

For the small angle MSW solution, a recent result for the solar neutrino requires |2j 
4c 2 „s 2 c 2 ~ (0.3 - 1) x 10~ 2 , which leads to 

Sl ~ 0.04 -0.07 (22) 

We shall next discuss the masses of the neutrinos. In general, the loop-level mass is 
suppressed compared to the tree-level mass simply due to the loop factors, giving a huge 
hierarchy among the neutrino masses. The atmospheric neutrino indicates the mass of the 
heaviest neutrino in the range of Eq. flT9|). The small angle MSW solution to the solar 
neutrino requires 

Am| MSW ~ (0.4 - 1) x 10" 5 eV 2 . (23) 

To explain the atmospheric neutrino anomaly and the solar neutrino problem simultaneously 
in this framework of the R-parity violation, one has to make the tree-level mass relatively 
small, which requires a partial cancellation in the AS 3 tan/3 + Am^ term. 

In the following, we shall demonstrate this fine tuning in the context of the GMSB [ IE, IS 



Similar conclusions will be obtained for other mechanisms of supersymmetry-breaking medi- 
ation, such as the gravity mediation. In a simple GMSB model where the messenger quark 
belongs to a pair of 5 and 5 representations, the masses of squarks and gauginos are char- 
acterized by one parameter, the effective SUSY breaking scale A. The universal trilinear A 
term as well as the universal B term are assumed to be generated. All these parameters 
are given at some energy scale, which we identify rather arbitrarily with A (= 10 5 GeV) in 
the following analysis. We also introduce the supersymmetric /i term and /Zj (i =1, 2, 3) as 
independent parameters. Given tan/?, defined as the ratio of the two VEVs of the neutral 
Higgs bosons, the conditions of the radiative electroweak symmetry breaking fix the values 
of \x and B. The value S3 is determined so that the heaviest neutrino mass lies in the region 
suggested by the atmospheric neutrino anomaly. Here we take its representative value to 

8 



be m V3 = 0.06 eV. As was explained earlier the other parameters S2 and si characterize the 
neutrino mixing angles. 

In Figs. 1 and 2, we plot contours of the ratio of the two neutrino masses, R = m U2 /m U3 , 
in the tan/? -A plane. Fig. 1 represents the case where /i < and Fig. 2 corresponds to the 
case of \x > 0. One finds that the preferred ratio R ~ 0.03 — 0.07 for the MSW solution to the 
solar neutrino is realized in small regions of the parameter space. We checked numerically 
that in these regions a cancellation among the two contributions to the VEV of the sneutrino 
indeed occurs. 

We also examined the value of S3. We found that the typical values of the S3 is large, 
compared to a naive expectation of (m^/m^) 1 / 2 ~ 10 -6 . Especially in the region where 
the MSW solution is obtained, we checked that the value of S3 can be ~ 10~ 3 — 10~ 5 . These 
observations confirm that the alignment of the vacua works at least partially so that one does 
not need to assume a priori a very small S3 to obtain small neutrino masses. Of course to 
simultaneously obtain the solutions of the atmospheric neutrino as well as the solar neutrino, 
we need the additional fine tuning to further suppress the sneutrino VEV. 

Let us now turn to the question of how one can test this scenario at future colliders. 
When the R-parity is violated, the lightest superparticle (LSP) which is also the lightest 
R-parity odd particle is no longer stable, but decays to ordinary particles. In this paper 
we want to discuss patterns of the decays of the LSP. Here we assume that the LSP is a 
neutralino.0 Indeed in the model described above a bino-dominant neutralino becomes the 
LSP. For simplicity we also assume that it is heavier than the W boson. The LSP then 
decays to the W (Z) boson and a charged lepton (a neutrino) through the sneutrino VEV. 
Or it decays to three body decays via scalar lepton exchanges using the trilinear R-parity 



2 In the GMSB models, actually the LSP is the gravitino. However, one finds that the decay of 
the neutralino into the gravitino is negligibly small, and thus one can concentrate on the R-parity 
violating decay of the LSP among the MSSM superparticles. 



violating couplings. In either case, we checked that the decay length is short enough for the 
LSP to decay inside a detector for typical parameters discussed earlier. In Fig. 3, we exhibit 
the ratio of the partial width of the two body decay which involves the vector boson to that 
of the three body decays, T(xi — > three leptons)/r(%5 — ► Wli,Zui). We find that for the 
MSW solution the two body decays are generally dominant decay modes, though the decays 
through the trilinear Yukawa couplings may not be negligibly small. 

Detailed analysis of the decays of the LSP may provide a test of the scenario we are 
discussing here. Let us illustrate this again when the LSP is bino-like. For the two body 
decay to a W boson and a charged lepton, the generation structure of the charged lepton 
directly reflects the mixing angles in the neutrino sector, and thus 



r(x?-A*wo __ \u 3 ,\ 



tan 2 6. (24) 



r(x? - rW) |f/3.| 2 

This will give information on the mixing angle of the atmospheric neutrino in this R-parity 
violation scenario. Note that the above relation is very generic, which is valid even when 
one relaxes the lepton-Higgs universality among the soft masses. 

On the other hand, the three body decays occur through the following Yukawa couplings, 

Y T s 2 s 3 ci(u T fi - u^t)t c + Y T s 2 s 3 si(u T e - u e r)r c (25) 

The neutralino LSP, which is in our case bino-like, couples most strongly to the right-handed 
stau and tau lepton and thus the decay occurs via the right-handed stau exchange and the 
final state will be either rrz/j or rv T li where the generation suffix i represents the first two 
generations. The latter final state is interesting because the tagging of the charged leptons 
in the first two generations will enable us to determine one parameter s\. In fact we find 

r(x?^rz/e) 2 

m^^fi 9l - {) 

This may provide us with an important information on our scenario where a large Si is 
excluded by the combined use of the CHOOZ experiment and the atmospheric neutrino. 

To summarize, we have discussed the case of bilinear R-parity violation as a source of 
neutrino masses and mixing. We found that the resulting mixing matrix of the neutrinos 

10 



has a very special pattern. This leads us to conclude that the large mixing angle solutions 
to the solar neutrino problem are ruled out when the CHOOZ result and the atmospheric 
neutrino data are combined together. Furthermore the relatively less hierarchical structure 
of the neutrino masses in this case are obtained if the soft SUSY breaking masses are suitably 
tuned to give small VEV for sneutrinos. Finally we argued that the analysis of the decays 
of the neutralino may reveal information on the mixing angles among the neutrinos. This 
scenario is testable at both neutrino oscillation experiments, e.g. SuperKamiokande pO 



SNO [pi]] , and KamLAND ||22| , and collider experiments in future. 
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FIG. 1. Contour plot of R = m U2 /m v . A on the trilinear coupling ^4(GeV)-tan (5 plane for [i < 0. 
The dashed line corresponds to R = 0.07 while the dot-dashed line corresponds to R = 0.03. 
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FIG. 2. The same figure as Fig. 1 except that \i > 0. 
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FIG. 3. The ratio of the decay widths T(xi ~* three leptons)/T(xi — ► Wk,Zui) versus the 
ratio of the neutrino masses m U2 /m U3 . The shaded region is obtained when R = 0.03 — 0.07 is 
imposed. 
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